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We reformulate boundary conditions for axisymmetric codimension-2 braneworlds in a way which 
is applicable to linear perturbation with various gauge conditions. Our interest is in the thin brane 
limit and thus this scheme assumes that the perturbations are also axisymmetric and that the 
surface energy-momentum tensor of the brane is proportional to its induced metric. An advantage 
of our scheme is that it allows much more freedom for convenient coordinate choices than the other 
^— ^ methods. This is because in our scheme, the coordinate system in the bulk and that on the brane 

■ are completely disentangled. Therefore, the latter does not need to be a subset of the former and 

^s^j ' the brane does not need to stay at a fixed bulk coordinate position. The boundary condition is 

manifestly doubly covariant: it is invariant under gauge transformations in the bulk and at the 
same time covariant under those on the brane. We take advantage of the double covariance when 
we analyze the linear perturbation of a particular model of six-dimensional braneworld with warped 
flux compactification. 
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INTRODUCTION 



O 

The idea of braneworld [1, 2, 3, 4, 5] has been attracting a great deal of physical interest as a scenario of our 
universe in higher-dimensional theories. In this scenario our universe is thought to be a timelike surface, called brane, 
embedded in a higher-dimensional bulk spacetime. The view of our universe as a 'membrane' has been inspired by 
the string theory and each braneworld scenario is considered to reflect some of essential features of the fundamental 
theory. In this sense, it is expected that various predictions from braneworld models should serve as tests to determine 
Oh! if the underlying theory is the true theory describing the universe or not. In particular, possible deviation of gravity 
from the four-dimensional theory has been considered as one of the most promising probes of the higher-dimensional 
theory. 

To compactify the extra dimensions is a key concept incorporated in higher-dimensional theories. In general, 
properties of gravity deviate from the prediction of the four-dimensional Einstein theory at scales shorter than the 
compactification scale. Experimental bound on the short distance deviation of gravity is not so strong and, in the 
context of braneworld scenarios, still allows the size of extra dimensions up to around 0.1 mm long [6, 7]. To be more 
precise, one should consider this as a bound on the mass of moduli, i.e., four-dimensional scalar fields representing 
deformations of extra dimensions. Even though the size of extra dimensions is smaller than 0.1 mm, light moduli 
would spoil the prediction of the four-dimensional Einstein theory at longer scales if the mass of the moduli were 
lighter than (0.1 mm) , For this reason, there have been tremendous efforts to stabilize all moduli in string theory 
[8] (for review, see [9] and references therein) and this has been one of the most important issues in string theory. 

As the opposite extreme, another key in the study of higher-dimensional theories is the expansion dynamics of the 
four-dimensional universe, in which the unusual behavior of gravity can in principle be observed. The non-standard 
expansion rate of the universe predicted by braneworld cosmologies serves as a test of the theories. Furthermore, one 
of the most fundamental problems assigned to the fundamental theory is realization of cosmological inflation [10, 11]. 
In this context, the use of the warped flux compactification has recently made it possible to investigate practical 
inflationary scenarios in the framework of the string theory [8, 12] (see [13] and references therein for review), and 
this can be considered as one of the most important achievements in very recent years. 
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For these reasons, the study of braneworlds has been getting more and more important. 

Now, the principle of general covariancc is one of the most essential concepts in modern physics: laws of physics 
must be independent of a priori geometry and the choice of coordinates, and, thus, covariant under general coordinate 
transformations. In the context of braneworld scenarios, there are indeed two types of coordinate transformation. 
One is on the brane and the other in the bulk. Laws of physics in our world on the brane must be doubly covariant, 
i.e., covariant under both types of coordinate transformation. To be more precise, they must be covariant under the 
coordinate transformation on the brane and invariant under the coordinate transformation in the bulk. 

Physically speaking, a brane representing our universe is not just a mathematical surface but must be manifested as 
an extended object or a soliton in higher-dimensional theories. Therefore it is in general necessary to specify properties 
(e.g., internal structure) of the extended object in order to describe dynamics of the braneworld. Fortunately, in the 
codimcnsion-1 case, it is possible to take the limit of zero thickness and obtain low-energy description of the dynamics 
without referring to microscopic properties. The thin limit can be taken also in the codimcnsion-2 case if the spacctime 
around a brane is axisymmetric and if the surface energy-momentum tensor includes vacuum energy, or tension, only. 
With codimension more than 2, it is not possible to take the thin limit, keeping brane tension and gravitational 
coupling constant finite. Indeed, for codimension more than 2, if we took the limit of zero thickness then gravity 
around the object would become so strong that it would form a black-hole horizon or a singularity. Therefore, for 
the description of a brane with codimension more than 2, one must keep the brane thickness finite and specify the 
internal structure within the thickness in order to describe the dynamics of the brane. From now on we shall restrict 
our attention to situations in which we can take the thin limit: the general cases with codimension 1 and the special 
case with codimension 2. 

In the codimension-1 case, Israel's junction condition [14] describes low-energy dynamics of braneworlds. It is beau- 
tifully formulated in terms of geometrical quantities (induced metric and extrinsic curvature) and matter quantities 
(surface energy-momentum tensor) on the hypersurface. However, a general scheme of perturbation of junction con- 
dition had not been formulated in a doubly covariant fashion until [15]. One of the obstacles was, to our knowledge, 
the fact that the coordinate position of the hypersurface changes under coordinate transformation in the bulk. For 
this reason, it had been common to use a part of gauge freedom to fix the coordinate position of the brane. This 
certainly simplifies the analysis in some cases, but it is a pity that the double covariance is lost. In [15] a doubly 
covariant perturbation scheme was given for the junction condition and it was applied to cosmological perturbation 
in the Randall-Sundrum braneworld. 

The purpose of this paper is to extend the doubly covariant perturbation scheme of Refs. [15, 16] to the codimension- 
2 case. As already stated, we shall restrict our attention to the situation in which the bulk geometry around the brane 
is axisymmetric and the surface energy-momentum tensor includes vacuum energy only. An advantage of the doubly 
covariant formulation is that we have more freedom for the choice of coordinates in the bulk than other formulations 
without double covariance. We shall take this advantage when we analyze the boundary condition for perturbation 
of six-dimensional warped flux compactification. 

As an application of our formalism, we investigate a six-dimensional braneworld model. In this model one or 
two 3-branes with tension are embedded in an axisymmetric warped spacetime, and the volume of compactified two 
dimensions is stabilized by virtue of a magnetic flux. The world- volume of 3-branes are 4-dimensional timelike surfaces, 
and one of them is considered as our universe. Specifically, we consider geometry of the form 

dr 2 

g MN dx M dx N = r 2 q„„{x)dx»dx v + — + f(r)dcf> 2 , 
A M dx M = A(r)d<(>, 

where q^ v is a four-dimensional metric 1 . A brane is located at r — rt>, where r^ is a zero of f(r), and thus q^ v 
actually represents the induced metric on the brane up to the constant conformal factor r\. The region with f(r) > 
represents the physical domain of the extra dimensions, while the region with /(r) < is unphysical and forbidden. 

This model has many interesting features. Indeed, within the ansatz (1.1), we can find a family of exact solutions 
of a braneworld with stabilized extra-dimensions. As already mentioned, one of the most important issues in higher- 
dimensional theories is moduli stabilization. For instance, in the Randall-Sundrum type-1 model, the modulus 
associated with the distance between the two branes acts as an additional degree of freedom appearing in low-energy 
description of gravity on the branes. As a result, the ordinary four-dimensional Einstein gravity would not recover, 
unless this modulus, called radion, is stabilized by some non-trivial mechanisms such as a bulk scalar field with suitable 
potentials. Unfortunately, no simple exact solution with a stabilized extra-dimension has been known in the context of 
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Randall-Sundrum type, five-dimensional braneworlds. On the other hand, with two (or more) extra-dimensions, one 
can introduce a U(l) gauge field (or an anti-symmetric field) to help stabilize the volume of extra dimensions. This 
is because all indices of non-vanishing components of the field strength can be entirely within the extra dimensions. 
The family of exact solutions is given by [17, 18, 19, 20] 

A-6 _ 2 Mb 



1 ~ ~ 10 ~ 73 1276' 2 ^ 

A = 373' 

and q^ v is set to be the metric of a four-dimensional Einstein space with the four-dimensional effective cosmological 
constant 3. Here, A 6 > is the bulk cosmological constant, fib < is a integration constant which describes the 
strength of the bulk Weyl tensor, and b corresponds to the magnetic charge of the U(l) field. The induced metric 
r^q^u on the brane is also an Einstein space but with the effective four-dimensional cosmological constant 3/r 2 . In 
[20] , the authors considered this family of solutions with q^ being the de Sitter metric and showed that the behavior 
of gravity on the brane recovers four-dimensional one at low energies. 

It is also possible to consider situations in which the effective cosmological constant of the induced metric on the 
brane vanishes. By taking the limit rb — > 00 and properly rescaling physical quantities, the above solution becomes 



t A 6 2 Hb 



A ~ 3r3' 



10 ' r 3 12r 6 

b 



(1.3) 



and q^ v is the metric of a four-dimensional Einstein space with a vanishing effective cosmological constant. In [21], it 
was shown that, at least when q^ is the Minkowski metric, this solution is stable against linear perturbations. Thus, 
the four-dimensional Einstein theory should be recovered in the linear and sufficiently low-energy regime. In that 
paper, we also showed relevant Kaluza-Klein (KK) spectra of the perturbations 2 . 

It is worth mentioning that these families of solutions exhibit some essential features of the warped flux compact- 
ification, including (i) warped extra-dimensions, (ii) moduli stabilization by a magnetic flux, and (iii) existence of 
brane(s). The warped flux compactification has been playing important roles in the studies of moduli stabilization 
in the framework of string theories. Therefore it is worthwhile to investigate these families of solutions to better 
understand general properties of gravity in warped flux compactification. 

There is also a variety of six-dimensional braneworld models which have similar structures to ours. Readers who 
are interested in such models should refer to, e.g., [22] and references therein. For other recent progresses, see also 
[23, 24, 25, 26, 27]. 

The rest of this paper is organized as follows. In section II, the formalism of the boundary condition for codimension- 
2 braneworlds is given. The boundary condition is doubly covariant with respect to coordinate transformations in the 
bulk and on the brane. In section III, application of our formalism to the above six-dimensional model is presented. 
In particular, it is analytically shown that in the limit where the warp factor is unity, the braneworld model is stable 
against linear perturbations. We summarize the result and discuss our future works in section IV. 



II. FORMALISM 



As already explained, the main subject of this paper is the dynamics of a brane with codimension 2. We would like 
to extract universal low-energy physics as much as possible without specifying microphysical properties associated 
with the internal structure of the brane. For this purpose we would like to take the thin limit, i.e., the limit in 
which the thickness of the brane is sufficiently small compared with physical scales of interest. This limit can be 
considered as just a low-energy limit, but it is known that rigorous treatment of this limit does not allow us to analyze 
a general setup but requires us to consider rather restricted situations. In particular, we have to assume that there is 
an axisymmetry in the extra dimensions at least in a neighborhood of the world-volume of the brane, that the brane is 
at the center of the symmetry and that the surface energy-momentum tensor includes tension only, i.e., a component 



2 In [21] it was concluded that there are zero modes in the tensor sector only. Actually, there are zero modes also in the vector sector as 
shown in section IIIB2 of the present paper. The existence of the vector zero modes does not change the main conclusion of [21] that 
the background is stable against linear perturbations and that the four-dimensional Einstein theory should be recovered at low energy. 
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proportional to the induced metric. Regardless of these restrictions, the thin limit has been useful for the analysis of 
codimcnsion-2 objects such as strings in four-dimensions and branes in higher dimensions. 

In this paper we shall still take these assumptions since we are interested in low-energy physics independent of 
the internal structure of the brane. The boundary condition at the brane is conveniently described by the following 
well-known relation between the deficit angle 5 and the brane tension a: 

' ° (2.1) 



MS+2 



where Mu + 2 is the Planck scale in the bulk spacetime. The deficit angle 5 in the left hand side of this formula can 
be easily calculated if the equation determining the brane position is explicitly written as w — Wb, where w is a bulk 
coordinate perpendicular to the world-volume of the brane and % is a constant. 

However, it is not always convenient to use such a coordinate system that one of the bulk coordinates is perpendicular 
to the world-volume of the brane and that the brane position is specified by a fixed value of this radial coordinate. 
Actually, in a coordinate system natural and/or convenient for the analysis of the bulk geometry, the brane position is 
often specified not by a fixed coordinate value but by parametric equations. In some cases including the examples in 
section III, the parametric equations even involve unknown functions, which arc determined only after solving some 
dynamical equations. Therefore it would be more convenient if there were a simple way to calculate the deficit angle 
in general bulk coordinates. 

In the doubly-covariant formalism developed in this section, the deficit angle can be calculated directly in any 
coordinate system and, thus, it allows more freedom for the choice of bulk coordinates. We shall take the full 
advantage of this freedom when we seek the boundary condition for perturbations of six-dimensional brane-worlds in 
section III. 



A. Deficit angle and boundary condition for axisymmetric configurations 

Let us consider a (D + 2)-dimensional bulk spacetime (M,g) and a D-dimensional world- volume B of a (D — l)-brane 
embedded in it 3 . The embedding is specified by the set of functions Z = {Z M (y)} as 

x M = Z M (y), (2.2) 

where x — {x M } denotes coordinates in the bulk and y — {y^} denotes coordinates on the world- volume of the brane. 
Hereafter, capital Latin indices run over all the dimensions M, N, . . . — 0, . . . , D + 1, and Greek and small Latin 
indices run, respectively, over the brane and the extra dimensions, i.e., ^,v, . . . = 0, ...,£)— 1 and a, b, . . . = D, D + 1. 

A neighborhood of B can be foliated by a 2-parameter family of D-dimensional timelike surfaces. Let us denote the 
corresponding embeddings by x M = Z^ (y), where £ = {Ci j C2} parameterizes the D-dimensional timelike surfaces, 
and suppose that Z^ =0 M (y) = Z M (y). The foliation is necessary just to make the geometrical quantities such as e p M 
and q^v (see below for their definitions) well-defined off the brane. The final expression of the boundary condition 
does not depend on the way the neighborhood of B is foliated, and is written in terms of local quantities on the brane 
only. 

Note that the coordinate system on the world- volume of the brane does not need to be a subset of that in the bulk. 
Thus, the coordinate transformations 

X M ^x' M =F M {x\ y^y>» = f»(y) (2.3) 

in the bulk and on the world- volume of the brane, respectively, are disentangled and are independently defined. Since 
the physical position of the brane is independent of the choice of coordinates, the embedding functions transform 
respectively as 

Z ( M (y) Z< M (y) = F M (Z ( (yj), Z ( M (y) Z<> M (y) = Z c {r\y)) (2.4) 

under these transformations, where / _1 is the inverse of the map / : y^ 1— > f^{y). Because of this property, the 
tangent vectors e p M defined by 



A p-brane is a p-dimcnsional extended object and, thus, its world-volume is a (p + l)-dimensional timelike surface. 
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transform as vectors under the bulk coordinate transformation. As a consequence, the product e^ M e u N gMN trans- 
forms as a scalar under the bulk coordinate transformation and, thus, the induced metric 

% v = en M ev N 9MN\ x=z{y) (2.6) 

is actually invariant under the bulk coordinate transformation. On the other hand, under the coordinate transforma- 
tion on the world-volume of the brane, the induced metric g M „ does indeed transform as a tensor. 

Now we assume that the two-dimensional extra space is axisymmetric and introduce the corresponding spacelike 
Killing vector field (p M in the (D + 2)-dimensional bulk spacetime and angular coordinate <j>. They satisfy ip M — 
(d/d(j)) M and Vm<Pn + VjvVm = 0, where Vm is the covariant derivative associated with gMN- We suppose that a 
brane is placed at the axis of the axisymmetry. Namely, we impose the condition 

L\ x =z (v) = 0, (2-7) 

where 



9 ± L \x=z {v) = -mr, (a-") 



L = y/g MN ip M ip N . (2.8) 

In the beginning of this subsection we introduced a 2-parameter family of D-dimcnsional timelike surfaces x M — 
Z^ M (y) to foliate a neighborhood of the world- volume B of the brane. We suppose that this also respects the 
axisymmetry, i.e., each surface in the family is orthogonal to ip: 

g M Ne^ M ip N = 0. (2.9) 

We shall relate the deficit angle to the normal derivative of the scalar L. For this purpose we introduce a unit 
vector field n M normal to the surface x M = Z^ M (y) and the Killing vector ip M : 

n M y M = n M e^ M = 0, 
riMn M = 1. 

These conditions determine n M up to the overall sign, which is fixed by demanding that n M should direct from 
the brane to the bulk. Note that n M is well-defined not only on the world-volume B of the brane but also in a 
neighborhood of it. The deficit angle 5 can now be written in a doubly-covariant way in terms of the embedding 
functions Z M (y), the unit norm n M and the scalar L: 

2tt - 5 

where d± = n M Om is the normal derivative and Atp is the period of the angular coordinate <f>, i.e., <f> ~ <f> + A<^. 

In summary, the boundary condition at the brane is given by (2.7) and (2.1), where the deficit angle 5 is given by 
the doubly-covariant formula (2.11). 

The doubly covariant boundary condition written above in terms of L is easy to calculate. Actually, we show an 
explicit result for a concrete example of background 6D warped flux compactification in section III. However, when 
we perform perturbative analysis around some background and seek the corresponding perturbative expansion of the 
boundary condition, this form is not most convenient since the definition (2.8) of L involves a square-root and thus 
is not an analytic function of metric components. Therefore, for the purpose of perturbative analysis, it is more 
convenient to use the following alternative expression: 

[y) = 0, (2.12) 
lz(„) - 0. ( 2 - 13 ) 

where S is related to the brane tension a by (2.1). It is evident that £0,1,2 are analytic functions of the metric 
components. Therefore, it is easy to perform perturbative expansion. 

One must, however, be aware of the fact that it is only up to overall sign oid±L\ x — z / y \ that this form of the boundary 
condition is equivalent to the original boundary condition. Actually this does not matter if £ 2 does not vanish for the 
background and if perturbation is small enough in the sense that the background value of 2n — 5 = 2n — cr/M£ +2 and 
the corresponding perturbed value have the same sign. This is the case unless the background value a is fine-tuned to 
the special value 2nM^ +2 with an extremely high accuracy. Anyway, if a had such a large value (<~ 2ttM£ +2 ) then 
the effective low-energy description would break down. Thus, while we can apply the original form of the boundary 
condition ((2.7) and (2.11) with (2.1)) to the background, we can use the alternative form ((2.12), (2.13) and (2.14) 
with (2.1)) for perturbation. 





L 2 \ 


1 X=. 






c 2 = 


dlL 2 \ 
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B. Linear perturbation 

The formalism developed in the previous subsection displays its real ability when used in perturbative analysis. 
Since the doubly-covariant boundary condition is written in terms of scalar quantities which are analytic functions of 
metric components, it is manifest how to expand it with respect to the perturbation. In this subsection we consider 
metric perturbation and perturbation of the embedding functions (representing the position of the brane), and expand 
the boundary condition up to the linearized order. 

We still assume axisymmetry and, thus, we consider axisymmetric perturbations only. 



1. Perturbed boundary condition 

Here we consider linearly perturbed geometry. As has been stressed, our purpose is to formulate everything in a 
manner which is covariant both in the bulk and on the brane. Since the coordinate transformation in the bulk changes 
the coordinate position of the brane, the double covariance requires us to consider perturbation of the brane position 
as well as metric perturbation. The most general perturbations of the background geometry and the position of the 
brane are 

9MN=g {0) MN + S 9 MN, V M = ^ M +6tp M , 

Z M (y) = Z^ M (y) + SZ M (y), 

where, throughout this paper, background quantities are recognized by the script '(0)'. Other linearly perturbed 
quantities are calculated in terms of SgMN, S(p M and 5Z M , and indices of perturbed tensors are raised and lowered 
by g^ MN and g^MN, respectively. Note, for example, that g MN = g^ MN — Sg MN up to the linear order. 

In the following we shall express the linear perturbations of £0,1.2 in terms of SgMN, S(p M and SZ M , where £0,1.2 
are defined by (2.12), (2.13), and (2.14). 

Before proceeding, we here make a comment on the perturbed Killing vector ip M . It is always possible to set 
Sip M — by using the gauge freedom and to set S(A(p) = by rescaling the angular coordinate (p. In the following 
analysis, however, while we shall set 5(A<j>) = 0, we shall suppose that dip M is not necessarily vanishing. This is 
because, in order to explicitly confirm that our perturbative scheme maintains the double covariance, we need to 
include 5ip M and its gauge transformation. Thus, we keep Sip M until section II B 3, where the perturbative double 
covariance is confirmed, but always set S(A</>) — 0. The perturbation Sip M of the Killing vector should satisfy the 
linearized Killing equation 

= 5 {Vm^n + Vjv^m) 

= -25T P MN ^ P + V<V M 6<p N + V (0) jv^m + V (0 W (0)P <S 9 PAr + V (0 W (0)P <W, (2.16) 

and should be orthogonal to the surface x M = Z M (y) so that 

= <5 (Sg MN ^ M e, N ) = 6g MN <pM M ef N + 6<p M ef M + <p™ M 6e™ , (2.17) 

where ST^ N = g(°) p Q(y(°) jv^fl'QM + V^M^fl'jVQ — V^°^q^mjv)/2 is the linear perturbation of the Christoffel symbol, 
and Se^ is that of the tangent vector: 

e M = e( Z)M + 8e M_ (218) 

(We shall give a concrete expression for 5e^ M in (2.24).) Needless to say, when we apply our formalism to concrete 
problems, it is more convenient to set 8ip M = by a gauge transformation than keeping it non- vanishing throughout. 
Therefore, in section II B 4 we shall simplify the perturbed boundary condition by setting Sip M = for later conve- 
nience. In section III, when we analyze perturbations of a particular six-dimensional braneworld model, we shall use 
the simplified expression. 

Let us now go back to the formulation of the perturbative scheme. 

First, we calculate the linear perturbation of the induced metric. The perturbation Sq^ is defined in the following 
manifest way: 

q lt »{v)=q$(v) + 8q liV (y), (2.19) 

where, up to the linear order, 

^{y) = 9MN e,V\ x=z{y) = 9M N e,V\ + SZ^d L (g^ MN e^e^) \ (2.20) 
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and 



JO) 



(0)MJO)N 



x=Z(°)(y) 



(2.21) 



Here, in order to get the last expression we have used the fact that gMNe^ M ' e v N is a scalar in the bulk and kept only 
terms up to the linear order. Thus we obtain 



SqM = S (9M N e, M e v N ) + 5Z L d L U a \ IN ef M e^ N ) 

x=ZW(y) \ / 

In general, linear perturbation of a quantity of the form X\ x=Z ( y ) is evaluated as 



x=Z(y) 



SX 



x=Z(°)(y) 



+ sz L d L x^ 



x = Z(0)(y) 



(2.22) 



(2.23) 



where X is a scalar in the bulk. 

In order to evaluate the first term in the right hand side of (2.22), we need to express <5e M M . Noting that e M M and 



Mm 



in equation (2.18) should be at a common point (say, at x = Z or x = Z^), it is understood that the linear 



perturbation Se^ 1 can be written in terms of the Lie derivative as 

5e„ M = £ em SZ M = -£ sz e^ M , 



(2.24) 



where £ denotes the Lie derivative in the bulk. (See equation (2.5) for the definition of e M M .) This expression allows 
us to calculate the linear perturbation of the induced metric. Combining (2.24) and (2.22), we can calculate Sq l _ l „(y) 
as 



Sq„(y) = [(Sg MN e^ M e^ N + g^M N Se, M e^ N + MN ef M ' 5e u N ) + SZ L d L (g^ MN ef M e^ N ) 
= (Sg M N + £6zg {0) MN)eW M e£ 



x=Z(°1(y) 



x=Z(°)(y) 

(2.25) 



Second, let us calculate the linear perturbation of Cq defined in (2.12). By the aid of the general formula (2.23) 
with X = gMNf M f N , this is an easy task. 



SC 



o(y) = s ( 



M Nl 



x=Z(y) 



3 (0 W (0) V 0)Ar 



s ( 9M n<p m <p n ) + sz L d L (. 

2<pW M 5<p M + ^ M ^ N Sg MN + SZ L d L (V°> V 0) p) 



x=Z(°)(y) 



x=Z(o)(y) 



(2.26) 



Third, let us calculate the linear perturbation of C\ defined in (2.13). By applying the formula (2.23) to X = 
d±(gMN^p M ^P N ), we obtain 



S&iy) = 5 (c*_l (g M N<p M <p N )\ x=z{y) ) 

5{d^{gMN V M V N ))+5Z L d L [df ] (. 



3 (0 W (0) V 0)Ar 



)) 



x=ZW(y) 



(2.27) 



To proceed further, we need to calculate the perturbation of the normal derivative d± = g MN riMdN- The perturbation 
of riM is actually determined by the following orthonormality condition: 



= S (n M <p M ) = <5rW° )M + n<°> M 6<p M , 

= 6 (n M e, M ) = Sn M e^ M - n^ M £ sz e^ M , 

- S (n M n M ) = 2Sn M n^ M - n^ M n^ N Sg M N, 



(2.28) 



where we have defined Sum as the linear perturbation of tim with a lower index: % = m + Sum- From these 
conditions, Sum is uniquely decomposed as 



Sn M = U^n^SgpQn^M - n ™? 6 £ ^°> M + ^ s£ S z^ N ef u . 



i P (o)Q (p (o) Q 



(2.29) 
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Remember again that this is the perturbation of Um with the lower index and that just raising the index of this 
expression by </°) M M does not give the correct formula for the perturbation of n M . Indeed, the perturbation of n M 
is given by S(n M ) = Sn M — 5g MN n^> n and, thus, 



S(d ± )=(8n M -6g MN nW N )d M . 
Using this expression for S(d±_), equation (2.27) gives the following expression for 5C\\ 

5Cr{y) = [(Sn M - n^ N Sg MN ) d M (V°> V 0) p) + n^ M d M (2<p™ P &<p p + ^ (0) V° )Q ^pq) 



(2.30) 



(2.31) 



+SZ L d L (n^ M d M (pW V 0) p 
Fourth, perturbation of C2 defined in (2.14) is also calculated in a similar way as 
SC 2 (y) = S (dlL\ =z{y) ) 

(5n M - nW N Sg MN ) d M (n<® L d L (V 0) V 0) p)) +n^ L d L ((<5n M - n^ N 6g MN ) d M (V 0) V 0) p) ) 



(> )M d M ) 2 (2^(°) P ^ P + ^ >V° )Q <5.9pq) +5Z L d L ^ M d M ) 2 (V 0) V 0) p)) 



X=Z(°) (y) 



Finally, the perturbed boundary condition is obtained as 

SC = Sd = 5£ 2 = 0, 
where #£0,1,2 are given by (2.26), (2.31), and (2.32), respectively. 



(2.32) 
(2.33) 



2. Boundary condition for U(l) gauge potential 

If a U(l) gauge potential Am^x m is present in the bulk, one must impose a separate boundary condition on the 
U(l) field. Here the perturbed Z7 (1) gauge potential is denoted as 



A 



M 



A (0) 



M 



5A 



ii/- 



(2.34) 



When we analyze physical perturbation of a background configuration, one should keep the total magnetic flux 
unchanged since it is a conserved quantity. With the total magnetic flux being the same as the background value, 
one can set the perturbation of A = AM ( p M \ x =z( y ) to zero on all the boundaries for the following reasons. First, the 
y-dependence of A is forbidden by the regularity of F MN Fmn- Second, the y- independent perturbation of A is not 
forbidden by the regularity but can be gauged away on all boundaries simultaneously since the total magnetic flux is 
assumed to be the same as the background. Therefore, we can set the perturbation of A to be zero on all boundaries. 

To calculate the perturbation of this quantity, one must take into account the fact that the position Z of the brane 
is also perturbed in general. Thus, the boundary condition for the perturbation of the U(l ) field is 



5A = 0, 



(2.35) 



where 



5A{y)^5(A M ^\ x=z{y) 

= 's(A MV M )+SZ L d L (A^ M ^ M 



X = Z(0) (y) 



6A M ^ M + A< W + 6Z L d L (>W(°) M ) 



x=Z(°)(y) 



(2.36) 



Here, we have used the formula (2.23) with X = Am<£ to get the second line. 

Note that, in general, it is not possible to set Am<P M itself to zero on more than one boundaries simultaneously. 
Indeed, with a non-zero total magnetic flux, a gauge in which Am¥> M is zero on a boundary is not compatible with 
another gauge in which Am^P M is zero on another boundary. The total magnetic flux, which is gauge-invariant, can 
indeed be expressed in terms of the difference between values of Am^P M at different boundaries; it is not the value of 
Am<P M \ x =z(y) itself but the difference 5 A that can be set to zero on all branes simultaneously. 
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3. Double covariance 



The general formalism developed in section II A is doubly covariant in the sense that it is covariant under coordinate 
transformation on the branc and invariant under coordinate transformation in the bulk. Therefore, the corresponding 
boundary conditions for linear perturbations should have double covariance up to the linearized order. In the following 
we shall confirm this explicitly as a consistency check. 

What we would like to show as a consistency check is that the perturbed boundary condition is doubly covariant 
under infinitesimal coordinate transformations in the bulk and on the brane, where the two kinds of transformation 
are denoted as 



.m 



respectively. 

Perturbations in the bulk, (2.15) and (2.34), transform under (2.37) as 



(2.37) 
(2.38) 



SgMN 
Sp> M 

SA M 
SZ M 



SgMN — £$g^MN, 
S^ M - £^ M , 
5 A M - £ t A<® M , 
5Z M +£ M . 



(2.39) 



Thus, it immediately follows from the formula (2.26) that 



SC 



2^°) M (Sp> M - £^ M ) + V° )7V (<W - £,9 

5g MN + SZ L d L (p^ p ^ P ) 



MN 



(sz L + e L ) d L ( 



x=Z(°~>(y) 



2^°> M V Vf + ^ 0) V 0) ' 
SCo. 



X=Z<.0)(y) 



(2.40) 



This says that the boundary condition SCo = is invariant under the gauge transformation in the bulk (2.37). 
Similarly, we can show that the other boundary conditions SC1.2 = and 5 A = are invariant under (2.37). 
On the other hand, under the gauge transformation on the brane (2.38), 



fidMN\ x =Z(°)(y) — * SgMN\ x= Z(«>(y), 



M 



\x=Z(°)(y) 



Sip 



5A M \ x =z(°)(y) 
7M 



Sip 
SA 



Mi 



|x=Z(0)(j/)j 



M| x= z(°)(j/)5 
M 



(2.41) 



5Z M -» SZ M - CM, 0)M . 

Then, from the formula (2.25), we see that the perturbation of the induced metric on the brane transforms as 

Sq, v - (Sg MN + (£sz £ c ^) 9 {0) mn) ef M e^ N 
= (5g MN + £ sz gW MN ) e^ M e£ 



,(0)N 



x = Z(0)(y) 



x=z(«H v ) 
£?a {0) 



(2.42) 



where £ D is the D-dimcnsional Lie derivative defined on the brane. Similarly, it is shown that #£0,1,2 and SA 
transform under (2.38) as 



(2.43) 



SCo- 


^SC Q 


f D r (0) 
— L, Q 


Sd - 


■+ Sd 


- £fC^ 


SL 2 - 


-» sc 2 


P D r (0) 
— X C 2 


SA- 


■+ 5 A - 





Therefore, all equations in the perturbed boundary condition have the double covariance as they should. 
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Simplified expression 



We have developed a perturbation scheme for our doubly covariant boundary condition. It is given by 

££ = 5Ci = 5£ 2 =SA = 0, (2.44) 

where ££0,1,2 an d &4 are given by (2.26), (2.31), (2.32), and (2.36). As a consistency check, we have explicitly seen 
that they are indeed doubly covariant in the linearized level. 

So far, we have been keeping 5(p M since we had to take into account its gauge transformation in order to show 
the perturbative double covariance in section II B 3. On the other hand, when we apply our formalism to concrete 
problems, it is more convenient to set 5ip M = by a gauge transformation than keeping it non-vanishing throughout. 
With 5ip M = (and 5(A<f>) = 0), we have slightly simplified expressions for ££0,1.2 and 6A as follows. 



££0- 
££1 = 

££2 = 
SA = 



S 9MN + SZ L d L ^ p ^ P ) 



X =Z(0) (y) 



(5n M n^ N Sg MN ) d M (V°) V° } p) + n^ M d M (V°) V° )Q £<7pq) 



+5Z L 8 L (n^ M d M V° } p 



x=Z«»(y) 



(2.45) 



(2.46) 



5n M_ n (0) 



6g MN ) d M (n(° )L d L V 0) p)) + n^ L d L [(5n M n^ N Sg MN ) d M (, 



( 

SA M ^ M + 5Z L d L (A^ M ^ M ) 



y°>V 0) p 



)) 



=Z(o)(l,) 



X =Z(°) (y) 



(2.47) 
(2.48) 



The perturbation of the induced metric is given by (2.25). 

Note that the boundary conditions formulated in this section should be supplemented with regularity of physically 
relevant, geometrical quantities such as the Ricci scalar of the induced metric on the brane, the tetrad components 
of the bulk Weyl tensor evaluated on the brane, etc. This is because we are adopting the thin brane approximation 
and all we can and should trust is what is obtained within the validity of this approximation. If, e.g., the Ricci 
scalar of the induced metric were singular then our approximation would be invalidated. It is of course possible to 
regularize the singularity by introducing a finite thickness of the brane. However, in this case the natural cutoff of 
the low-energy effective theory is the inverse of the thickness, and in general the regularized 'would-be singularity' is 
not expected to be below the cutoff scale. This simply means that we need a UV completion, e.g., the microphysical 
description of the brane, to describe the physics of the regularized 'would-be singularity'. Therefore, in general we 
have two options: (i) to specify a fundamental theory such as string theory as a UV completion and go on, or (ii) 
to concentrate on modes which are within the validity of the effective theory. In the present approach we adopt the 
latter attitude, assuming the existence of a good UV completion but never using its properties. This is the reason 
why we adopt the thin brane approximation and require the regularities. 

In the next section, we shall apply these expressions to perturbations of the six-dimensional braneworld model 
explained in the introduction. 



III. APPLICATION TO CODIMENSION-2 BRANEWORLD 



Now wc apply the formalism developed in the previous section to a concrete example. We investigate perturbation 
of a six-dimensional braneworld model with warped flux compactification. Two families of background solutions were 
introduced in section I. The first family is given by (1.1) together with (1.2), where is a four-dimensional Einstein 
space with the effective cosmological constant 3. The second family is obtained as a limit of the first family and 
is specified by (1.1) with (1.3), where q^ v is an Einstein space with vanishing effective cosmological constant. As 
explained in the introduction, these families of solutions have many interesting aspects as a model of braneworld. 

For these backgrounds, the Killing vector is 

^ (0)M = <, (3.1) 

corresponding to the angular coordinate <f>. We denote the period of the angular coordinate by A<j>, i.e., <f> <~ <f) + A<f>. 
Since is written as 



= y/C = 



(3.2) 
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the boundary condition (2.7) for the background says that the unperturbed position r^ ' of the brane is a zero of the 
metric function f(r), as mentioned in the introduction. The unit normal to the brane is 



n (0)M = Tv/79 M 



(3.3) 



where the minus and plus signs are for the cases f'(r^) < and f'(r^) > 0, respectively 4 . Thus, wc obtain 

af LW= Ty /fd r y/f = T\f'(r). (3.4) 
The corresponding boundary condition (2.11) relates the parameters of the solution to the brane tension o\, as 



2tt - S h _ 
A<f> ' b ~M 4 ' 



(3.5) 



For appropriate values of brane tension, we can put one or two de Sitter or Minkowski branes at zeros of the function 
f(r) by imposing this boundary condition. For more details of the properties of the solutions, see [20]. 

Linear perturbation around these backgrounds is described by 5gMN, 5 Am and 8Z M . (We have set 5ip M = and 
S(Aip) — 0. See the third paragraph of section IIB 1, especially the sentences just after equation (2.18).) As already 
stated in the formulation, we still assume that the bulk spacetime has an axisymmetry and that the brane is at the 
center of the symmetry. Hence, SZ M can be written as 



SZ 



M 



5r b d, 



M 



(3.6) 



The boundary condition is given by (2.44), where 5 Co, 1,2 and 5A are calculated by the formulas (2.45), (2.46), (2.47). 
and (2.48) as 

SC = [6 g </,</, + f'Sr h ] r=r m , 
5d = 



-\f /2 f'5grr + y/fSffa + (v7/')'*rb 

-f 2 (Vff')'Sgrr - \ff (fSgrr)' + ^1 (v7<fo^) ' + (v 7 / (v7/')') ^ 



5C 2 

SA = [5A,p + A'5r h ] 



(3.7) 



Here, a prime denotes derivative with respect to r. The second-order derivative Sg'^ can actually be eliminated by 
using the equation of motion. 



A. Scalar perturbations in the 6D braneworld and the boundary condition 

The boundary condition given above is applicable to perturbations around both families of background solutions. 
In the following, we shall restrict our consideration to the second class of background solutions, i.e., (1.1) and (1.3), 
and set q^ v to be the Minkowski metric. Linear perturbations around background solutions in this sub-class were 
studied extensively in our previous work [21]. In particular, perturbative stability was shown and the Kaluza Klein 
spectra were obtained. 

In [21], the boundary condition for scalar- type perturbations was not derived from the first principle but obtained 
by assuming that a particular set of perturbation variables and their derivatives with respect to the coordinate r 
should be finite on the brane. However, it is not a priori clear which variables should be finite on the brane. Suppose 
we have two perturbation variables $i and $ 2 - It is also possible to use another combination, say r 2 $i and $2 + H&i- 
Indeed, there are in principle infinite possibilities for the choice of perturbation variables. The finiteness of which set 



4 Note that we have chosen the direction of n M so that it directs from the brane to the bulk. Thus, the sign convention here was chosen 
so that, when the bulk is bounded by two roots r± (r_ < r+) of /, the minus sign in (3.3) is for the larger root r+ and the plus sign 
is for the smaller root r_. We do not consider the degenerate case with /'(^b) = 0. Actually, within the family of solutions considered 
here, there is no regime of parameters that gives a degenerate root of / and compact extra-dimensions. 
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of perturbation variables should we impose? Even if we can somehow guess the answer to this question, there are in 
principle infinite possibilities for the choice of the coordinate with respect to which the derivatives of the perturbation 
variables are set to be finite on the brane. Therefore, it is crucial to justify the boundary condition used in the 
analysis of [21] from the first principle. This is what we shall do in the following. Moreover, in section IIIB we shall 
analytically show the stability of the background in a particular limit called the football-shape limit. 
Here we concentrate on scalar- type perturbations with the gauge choice adopted in [21]: 

dr 2 

g MN dx M dx N = (1 + ^Y)r 2 ^dx^dx" + 2h {L)(j> V {L)tl dx^d^ + [1 + ($ x + <t> 2 )Y}— + [1 - ($ x + 3$ 2 )F]/d0 2 , 



A M dx M = a r Ydr + (A + a^Y)d^>. 



(3.8) 



Here the metric components in the right-hand side have been already expanded by the harmonics of the Minkowski 
space summarized in Appendix A, and integration in the k space is understood although it is not written explicitly. 
For completeness, Appendix B includes detailed explanation of the gauge transformation leading to this gauge choice. 
The coefficients {$i, $2, ^(l)</>7 a r , a<f,} are functions of r and represent physical perturbations up to a residual gauge 
freedom C explained in Appendix B. We also expand 5r\, as 



I 



dkSr h Y. 



(3.9) 



Correspondingly, by using (3.7), we obtain the harmonic expansion of quantities 5£o,i,2 and 8A relevant for the 
perturbed boundary condition (2.44): 



SCo, 



dk5C 012 Y, SA 



dkSAY, 



(3.10) 



where 



SC = 
5C t = 

sc 2 = 

5A = 



-/($i+ 3$ 2 ) + f'Sr h 



„(t>) ' 



~V7/'(3*i + 7* 2 ) - f/ 2 (<S>[ + 3$ 2 ) + f2 ^ f " sr h 



(0) 



-if' 2 + 2//")(*i + 24> 2 ) - //'(3*i + 8$ 2 ) - / 2 ($'/ + 3$2) + (2/'/" + f f")5r h 
+ A'5r h 



(3.11) 



Thus, the perturbed boundary condition becomes 



££0,1,2 = SA = Q. 



(3.12) 



In the following we shall simplify the expression (3.11) and show that this boundary condition is actually equivalent 
to the naive boundary condition used in [21]. 

For this purpose, first, we shall use the field equations in the bulk. As shown in [21], we can algebraically solve a 
part of Einstein and Maxwell equations to express four of six variables in terms of the other two: 



A' 1 

* = $2, /i(L)0 = Cf, a r = — jh {h)cf> = -CA', a 4> = — 



(3.13) 



where C is a constant corresponding to the residual gauge represented by C in Appendix B. The remaining field 
equations reduce to 



4 77? 
*2 + -* , 2+2^7(*i+2* 2 )=0 ) 



f 5 
/ + r 



/ r 2 f 



(3.14) 
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where m 2 = —rj^k^ky. By using these equations of motion, <5£ 2 in (3.11) is reduced to 



SC 2 



.f 2 + / 4A, 



5m 
^2" 



+ 2/" $1 - 2.f' 2 + / 4A 



3m 2 



+ 4.f 



$ 2 



( /' 4/ (V - ^ ) <I> 2 ,2/'/" - //"»„ 



(3.15) 



Second, as explained in the end of section II, we require that geometrical quantities be regular on the branc. In 
particular, we require that linear perturbations of R, R mnp ®Rmnpq, V m RVmR, and V R R MNP ®\7 r R MNP q to 
be regular on the boundaries, where / vanishes. By using computer algebra packages such as GRTensor [28], we can 
express these geometrical quantities in terms of $i j2 and their derivatives. After eliminating second and higher order 
derivatives by (3.14), they are reduced to linear combinations of just <&i, 2 and $'i 2 , w i tn coefficients being functions of 
r. At this point, however, the matrix made of the coefficients does not have a regular inverse in the limit / — > 0. After 
some algebra, it is found that the linear transformation from 7^1,2.3,4 to linear perturbations of R, R MNP ® Rmnpq, 



V m RV m R, and V R R MNp Q\7 R R 



R-tlMNPQ 



is regular and has a regular inverse in the limit r 



„(°) 



where 7£i 2 3 4 are 



defined as Tlx = $2, ^2 = ^3 = /' $2 ~ (m 2 /2r 2 )^ 1 , and ft 4 = (/$i)'. Therefore, ^1,2,3,4 should be regular 
on the boundary. Since 7£i, 2,3,4 multiplied by positive powers of / vanish on the boundary, let us express <5£o,i,2 in 
terms of 7£i, 2,3,4 instead of $1,2 and 2 in order to identify terms remaining on the boundary. The result is 



SC a = 

sr.! = 

5C 2 = 



-TZ 2 - + f'Sr h 
3m 2 ^f f 



AO) ' 



2r 2 



/' ^)/v; 1/(2 ^) 1Z, - riff - /./"%»„ 



,<°> 



4A 6 + ^ + ^ +2 r-^>2-f 2 r + /f4A 



,-^+4/") )■/?: 



(3.16) 



Since 72-1,2,3,4 are regular on the brane, all terms in the above expressions are manifestly regular. Terms with positive 
powers of / all vanish on the brane and, thus, we obtain the following expressions after recovering $1 , 2 and $' L2 . 



S£ = 

sr.! = 

6C 2 = 



7*i + f'Sr h 



$1 



2V7 



„(0) ' 



Sr h 



4A 6 + ^ + ^ + 2/"j /'<!., - 2/'-<l» 2 - /'(./<I. l .)' + 2/'/"^V, ) 



(3.17) 



Actually, the four boundary conditions ££0,1,2 = SA = are not independent but only three of them are independent 
since is now written in terms of $i, 2 as in (3.13). Hence the independent conditions are 



/*l| r=r (0, - [2/'$ 2 + (/*!)'], 



,(0) 



0, 



(3.18) 



and 



5r b = 0. 



(3.19) 



The last equation means that the brane position in this gauge is not changed by the perturbation. Note, however, that 
we would not be able to know this fact unless we used the doubly covariant formalism developed in this paper. With 



the manifest replacements r 



(o) 



.(0) 



and <5rb — > Sr±, the boundary condition (3.18) is exactly what we obtained in 



[21] by naively setting coefficients of / _1 to be zero. On the other hand, in the present paper we have reached the 
same boundary condition from the first principle. Therefore, we have justified the treatment in [21]. 

For vector and tensor perturbations, the perturbed boundary condition (2.44) does not give any requirement since 
the quantities £0,1,2 and A are scalars. Nonetheless, as shown in [21], regularity of geometrical quantities on the brane 
gives sufficient boundary conditions. 
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B. Stability of the football-shape limit 

Next we consider a limit in which the shape of the extra dimension becomes locally a round 2-sphcrc. Taking the 
a = r_/r + — > 1 limit gives 

g (0) MNdx M dx N = rlrj^dx^dx" + R 2 fl ab dx a dx b , (3.20) 

where the constant ro represents the position of the two branes and R is the radius of the sphere related to the bulk 
cosmological constant as R = 1/V2A 6 . f2 a b is the metric of round sphere and <j> parameterizes its azimuthal angle 
whose period is 2n — S, i.e., the sphere has a deficit angle S = cr/M^, where a is the common tension of the branes. 
Although the coordinate distance r + — r_ between the branes vanishes and the bulk geometry seems to disappear, 
the proper distance remains finite. For the detail of this limit, see Appendix C. In the following we normalize the 
warp factor at the position of the branes and the radius of the sphere so that r = R = 1 . The geometry of the extra 
dimensions has a coordinate representation 

n ab dx a dx b ^^f- + j{w)d^ 1 

/H (3.21) 

A M dx M = A(w)dcp, 

in which the metric function and the gauge potential are specified as 

/ = l-w 2 , A = -w. (3.22) 

The background quantities are given as 

^(0)m =9 m L (0) = v /7-^ n (0)M = yT3^M ? (0)^ = 7W) 4°) =±L (3.23) 

Perturbations are defined on this background. Schematically, all the formalism developed above can be directly 
applied to this limit together with the following replacements 

(r,0) ->(«;,¥>), /->/, A^A. (3.24) 

In the following, we will show stability of all the types of perturbation for completeness. However, the boundary 
conditions obtained above are relevant for the scalar perturbations only. Thus if the reader is particularly interested 
in direct applications, it is recommended to skip to section IIIB3. As will be shown there, the boundary conditions 
for the scalar perturbations at the branes (3.18) are further reduced to 

[*i+2* 2 U ±1 =0. (3.25) 

1. Tensor perturbation 
The configuration with tensor perturbation is given by 

2 



g M Ndx M dx N = + fc ( T)T (T ) ^dx^dx" + + (1 - w 2 )d^ 2 , 

Audx M = —wdtp, 

where h(r) denotes physical perturbation. Tensor perturbation is gauge- invariant. 
The Einstein equations are reduced to 

[(1 - u; 2 )^)]' -iVtyT) = 0, (3.27) 

where the prime denotes derivative over w and the (dimensionlcss) Kaluza-Klein mass squared /i 2 is defined by 
/i 2 = —ri^k^ku. The solution is obtained as 

h {T) = cP v {w) + dQ v (w), (3.28) 

where P v and Q v are the Legendre functions of the first and second kind, respectively, and c and d are arbitrary 
constants. Their order v is expressed by above defined /j, 2 as 

„=Z±±VM±1. (3.29) 
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The condition imposed on the tensor perturbation is regularity on the boundaries w = ±1. By requiring the Ricci 
scalar of the induced metric and vielbein components of the Weyl tensor to be regular there, we find that fyx) and h'^ 
should be both regular on the boundaries. The set of independent regular quantities used here is listed in Appendix D. 
The regularity of /im at w = 1 requires d = 0, and that at w = — 1 does P v itself to be non-singular there, which 
is realized when v is non-negative integer: v = 0, 1,2, . . .. The mass spectrum of the tensor perturbation is then 
determined as 

fi 2 = is(is + l) [y = 0,1,2,...). (3.30) 

The spectrum of the first few KK modes is listed in Table I. Note that for simplicity we have been setting the radius 
of the extra dimensions as R = l/\/2k§ = 1, however, in general cases the mass spectrum is given by 2Aq^, 2 and this 
scaling applies also to the vector and scalar perturbations discussed in what follows. 



TABLE I: The KK spectrum of the tensor perturbation. 



Level v 





1 


2 


3 


4 


5 




Mass fi 2 





2 


6 


12 


20 


30 





There is a zero mode describing the four-dimensional gravitational wave and other modes have positive /z 2 hence 
we conclude that the spacetime is stable against linear tensor perturbation. 



(3.31) 



2. Vector perturbation 
The configuration with vector perturbation is given by 

g M Ndx M dx N = ^dx^dx" + 2h {T)a V {T) ^dx a + ^—^ + (1 - w 2 )^ 2 
A M dx M = a (T )V"(T) A1 da; A1 - wd(p, 

where {h^ w , ft.(x) v , O(t)} denote physical perturbations. Gauge fixing is summarized in Appendix B. 
First, the [iv components of the Einstein equation is solved to give 

k{T)w = 1 - w 2 ' < ^ 3 ' 32 ' ) 

where D is an arbitrary constant when /i 2 = —r/^k^k^ = 0, or D = when fj, 2 ^ 0. However, if D has a non-zero 
value, divergence of the metric components on the boundaries cannot be removed since no gauge freedom is remaining. 
Hence we conclude that D — 0. Remaining Einstein's equation and Maxwell equation reduce to 

(1 - w 2 )(h>( T)v - 2o' (T) ) + f fyi> = 0, 
[(1 - w 2 )a{ T) }' + h{ T)v + M 2 a ( T) = 0. 

By requiring the vielbein components of the ga uge field strength and those of the Weyl tensor to be regular on the 
boundaries w = ±1, we find that a(T)> Vl — w 2 a^ T ^ ft.(x) v /(l — w 2 ), and h'^ T ^ should be regular on the boundaries. 
The set of independent regular quantities used here is listed in Appendix D. 
In the case of n 2 = 0, h'^ T ^ is solved as 

h[ T)v = cP 1 (w) + dQ 1 (w), (3.34) 

where c and d are constants, and Pi and Qi are, respectively, Legendre polynomials of the first and second kind 
of order 1. Requiring the regularity of /i(x) ¥ ,/(l — w 2 ) and h'^ v at w = ±1, we obtain the zero mode of vector 
perturbation, 

h(T)<p — ci(l - w 2 ), a(T) = -c\w + c 2 , (3.35) 
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where c,'s are arbitrary constants corresponding to two degrees of freedom for the zero mode of vector perturbation. 
Next we consider Kaluza-Klein modes. When /i 2 7^ 0, introducing following variables makes the analysis easier: 

h W)<P = Tit 1 - w2 ) $ '(V)2, 0(T)=*(V)1- (3-36) 

Then the field equations (3.33) become 

[(1 - ™ 2 )$' (v)2 ]' + ^ 2 (<i>(v)2 - 2*(v)i) = 0, 
[(1 - ™ 2 )$' (v)1 ]' + (M 2 + 2)$(v)i - * ( v)2 = 0, 

where $(y)2 was remove d by the equations of motion. Eliminating $(v)i from these equations gives a single fourth- 
order differential equation for $(v)2 

[(1 - w 2 )^]' + \±E± = 0, (3.38) 

where 

^ ± = [(i-^ 2 )^ ( v )2 r+A T j> ( v )2; 

A± = ^ 2 + 1 ± + 1. 

The case of ,u 2 = —1/2 causes the eigenvalues A± degenerate and requires another treatment. If /x 2 7^ —1/2, the 
solution for $(y)2 is obtained as a linear combination of the Legendre functions and, from (3.37), $(v)i is readily 
solved. The general solutions are 

*(V)2 = c+P v+ + c-P v _ + d + Q v+ + d-Q v _, 

$ (v) i = [c+(A+ - M 2 )^ + + c_(A_ -f)P v _ +d+(\+ -m 2 )0, + +rf-(A_ -/i 2 )Q„_] , 

where c± , d-t are arbitrary constants and the values of order v± are determined by 

v±(v± + l) = X±. (3.41) 



Now $(v)2> ^(v)2' ^(v)ii an( i Vl — w2 ^[v)i are independent regular variables on the boundaries. At w = 1, 3>(v)2 is 
regular only if c? + + cL = whereas the regularity of $(v)i requires rf+(A + — fi 2 ) + gL(A_ — /1 2 ) = 0. Since A + 7^ A_ 
for /i 2 7^ —1/2, we obtain d± = 0. At io = —1, the regularity of the variables requires 

= c+ sin^+7r + c_ sini/_7T, 

, . x • (3-42) 

= c + A+ sin v + tt + c_ A_ sin V-it. 

With A + 7^ A_, this implies that non-trivial solutions can exist if v + e Z or e Z. We have to only consider 
non-negative i/±'s so that these conditions are explicitly written down as 



-1 + a/V + 5± V2/z 2 + l 
i/± = * = 0,1,2,.... (3.43) 

This condition gives the KK mass spectrum for the vector perturbation. The first few of them is listed in Table II. 



TABLE II: The KK spectrum of the vector perturbation. 



Level 


V- = 0, v+ = 1 


v+ = 2 


V- = 1 


v+ = 3 


v-=2 


v + = 4 




Mass \j? 





2(3 - y/3) 


4 


2(6 - V6) 


2(3 + y/3) 


2(10 - VW) 
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As mentioned earlier, the zero mode of vector perturbation has two degrees of freedom, thus there are v + = 1 and 
v_ = at /i 2 = 0. 

In the special case of yu 2 = —1/2 the eigenvalues A± degenerate and the order is given a value v= (VS- l)/2. The 
general solutions are given by 



*(V)2 = [cP v + cd v P v + dQ v + d'd u Q u } u=(x ^_ 1)/2 , 
*(v)i = [c^ + c'(V3 + <9„)P„ + dQ v + d'(VTi + d v )Q v 



(3.44) 



i/=(\/3-l)/2 



where c, c', d, d' are constants. At w = 1, they are regular only if d = d! = 0. Then, to remove divergence at w = — 1, 
we also obtain c = c' = 0. Therefore ^i 2 = —1/2 mode does not exist. 

All the KK modes have positive /z 2 hence we conclude the spacetime is stable against linear vector perturbation. 



3. Scalar perturbation 

Finally we investigate scalar perturbations. Now it is required to consider boundary conditions derived by our 
formalism. This in contrast with the other types of perturbations, for which the regularity of geometrical quantities 
was sufficient. 

The configuration with scalar perturbation is given by 

g M Ndx M dx N = (1 + VY^dafdaT + 2h {h)v V (L) ^d V + [1 + (*j + $ 2 )y] T _2 

+ [l-($ 1 +3$ 2 )y](l- W 2 )d^ 2 ( 3 - 45 ) 
AM<ix M = a w Ydw + (—w + a ip Y)d<p 

where we have already chosen an analogue of the longitudinal gauge so that {$i, $ 2 , ^, ^(L) VJ a w ,a v } denote physical 
perturbation. Note that there remains a residual gauge freedom in h^ v . Gauge fixing is summarized in Appendix B. 
From Einstein's equations, four of six variables are algebraically solved as 

* = 4>2, /*(L) V = C(l - w 2 ), a w = - X -h!'^ )v = C, a v = W ($ 1 +2$ 2 )-i(l- w 2 )cI>' 1 , (3.46) 

where C is a constant corresponding to the residual gauge, which is represented by C in Appendix B. Then the 
remaining field equations reduce to 

(l- u , 2 )$ 2 ' + ^($ 1 +2$ 2 )=0, (34?) 
(1 - w 2 )<5>'{ - 4w<f>[ - (2 - ^ 2 )$! - 2$ 2 = 0, 

where [i 2 = —rfk^kv. By requiring the Ricci scalar, Kretschmann scalar, vielbein components of the Weyl tensor to 
be regular on the boundaries w = ±1, we find that $ 2 , $ 2 , <&i, and \/i — w 2 $' x — (2w/Vl — w 2 )(<J>i + 2$ 2 ) should 
be regular on the boundaries. Multiplying the last one by \f \ — w 2 and given that both <i> 2 and $1 is regular at the 
boundaries, we deduce that (1 — w 2 )&[ is also regular. The set of independent regular quantities used here is listed 
in Appendix D. 

Following analysis needs the aid of the boundary condition (3.7). By the knowledge of the regular variables, we 
can drop terms consist of regular combination multiplied by positive powers of 1 — w 2 . Second order derivative of the 
metric component appearing in 8C2 is removed by the equation of motion (3.47). After such reductions, we obtain 
the boundary conditions as shown in (3.25), 

[*i + 2* 2 U ±1 =0, (3.48) 

and the regularity at w = ±1 yields a v \ w= ±i = 0. 

Now we show the absence of zero mode. When /j, 2 = 0, the bulk solution for (3.47) is 

$2 = ciw + c 2 , 

ciw 3 + 3c 2 w 2 + c 3 w + c 4 (3.49) 
1 3(1 -w 2 ) 
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where Cj's are arbitrary constants. This system has in fact no solution other than $1 = $2 = when the regularity 
at the boundaries and the boundary condition (3.48) are imposed. 

Next we consider Kaluza-Klein modes. Eliminating <I>i from the field equations (3.47) gives a single forth-order 
differential equation for $ 2 



where 



[(1 - w 2 )E' ± ]' + \ T E± = 0, 



E ± = [(1 - w 2 )& 2 ]' + A±<1> 2 , 



a± = [i 2 + 1 ± vV+T. 



(3.50) 
(3.51) 



The case of /x 2 = —1/3 causes the eigenvalues A± degenerate and requires another treatment. If [i 2 7^ —1/3, the 
solution for <J> 2 is obtained as a linear combination of the Legendre functions, and, from (3.47), $1 is readily solved. 
The general solutions are 



$2 = c + P v+ + c-P v _ + d + Q u+ + d-Q v _ , 
2 

$1 = — [c+(A+ -[i 2 - 2wd w )P v+ + c_(A_ - /1 2 - 2wd w )P v _ 
I- 1 

+d + (X + - /j, 2 - 2wd w )Q v+ + d_(A_ - [l 2 - 2wd w )Q v _], 
where c± , d± are arbitrary constants and the values of order v± are determined by 

v±{v± + 1) = A±. 



(3.52) 



(3.53) 



Again we insist that, for the scalar perturbation, it is not sufficient to impose regularity on the variables but the 
boundary condition is needed. The boundary condition for $1 and $2 is expanded around w = 1 as 



= [$1 + 2$ 2 ] w=1 
~ \(d+ + d_)(l - w)- 1 - -^[d+(\+ - 1) + d_(A_ - 1)] + 0(1 - w). 



(3.54) 



Since A + 7^ A_, this means d+ = gL = 0. Next we expand it around w = — 1, then 

0= [$i+2$ 2 ]„ = _ 1 

4 2 
~ — -(c + sin v + -k + c_ sin^_7r)(l + w)^ 1 H -[c + (A+ — 1) sinz/ + 7r + c_(A_ — 1) sinz/_7rl + 0(1 + w). (3.55) 

This, again with A+ 7^ A_, shows that non-trivial solutions can exist if v+ € Z or i/_ e Z. We have to only consider 
non-negative i/±'s so that these conditions are explicitly written down as 



-1 + \Up? + 5 ± 4^3m 2 + 1 
"± = V ~ 7: = 0,1,2,. 



(3.56) 



Substituting (3.52) into the expression for a v with a non-negative integer v = v±, we can explicitly examine a v \ w= ±x = 
0. Finally, the above condition for v gives the KK mass spectrum for scalar perturbation, which is listed in Table III. 



TABLE III: The KK spectrum of the scalar perturbation. 



Level 


v-=Q 


v+ = 2 


V- = l 


v+ = 3 


v-=2 


v+ = 4 




Mass \j? 


1 


(13 - v / 73)/2 


5 


(25 - v / 145)/2 


(13 + V73)/2 


(41 - v / 241)/2 





In the case of [i 2 = —1/3, the eigenvalues A± degenerate and the order is given a value v — (\/33 — 3)/6. The 
general solutions are given by 



$2 = [cP v + dd v P v + dQ v + d'd u Q u ] _ vm _ 3 , 



3c(l - 2wd w )P v + c' (\/33 + 3(1 - 2wd w )d^j P v 
+3d(l - 2u;0 1u )Q 1 , + d' (V33 + 3(1 - 2^)9,) Q„ 



(3.57) 



V33-3 : 

6 
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where c,c',d,d' are constants. At w = 1, they are regular only if d = d! = 0. Then, to remove divergence of at 
w = — 1, we also obtain c = c' = 0. Therefore /x 2 = —1/3 mode does not exist. 

All the KK modes have positive /j, 2 hence we conclude the spacetime is stable against linear scalar perturbation. 



In this paper, we have formulated an adequate boundary condition for linear perturbations of axisymmetric 
codimension-2 braneworlds, extending the doubly covariant formulation of [15] for codimension-1 braneworlds. An 
advantage of this formulation is that it allows much more freedom for convenient gauge choices than others. The 
boundary condition has been expanded in perturbations and resulted in equation (2.44), where &Co,i,2 and 8 A are 
given by (2.45), (2.46), (2.47), and (2.48). 

As an application, we have investigated linear perturbations of a six-dimensional braneworld model in section III. 
We have derived the boundary condition for scalar- type perturbations that was used in our previous analysis [21], 
where the dynamical stability of the model with Minkowski branes was shown numerically. In the present paper, 
moreover, we have analytically shown the stability of the system in the limit where the warp factor becomes unity. 

The present status of our understanding of gravity in codimension-2 braneworlds is still premature and indicates 
that more investigations are necessary. Our future work in this direction has at least two branches. One is the 
extension of our formalism to more general setup where the geometry in the bulk is not axisymmetric and/or where 
general matter contents other than tension are put on the branes. In particular, finite thickness of the brane will take 
important roles in this context [29]. It is expected that the extension will make it possible to reveal more properties 
of gravity in the warped flux compactification. The other is to investigate the nature of the geometry in the high- 
energy regime. The importance of the solutions with de Sitter branes comes from its connections to the inflationary 
early universe and/or the present accelerating expansion. In [20], a sign of instability was noticed when the Hubble 
expansion rate on the brane is sufficiently large. In a forthcoming publication we hope to investigate the possibility 
of such an instability and to understand its origin. It is also worthwhile to consider a black hole on the brane [30] 
and investigate its stability. 

Further developments are needed also in the context of quantum dynamics of the braneworld. In particular, 
quantization of matter contents is to give important predictions of inflation in the early universe, and the most direct 
way of quantization requires the action of the system. In this respect, it is preferable if we can obtain a doubly 
covariant action including not only fields on the brane and those in the bulk but also the brane position as an 
independent variable. It had actually been believed for a long time that the action principle of singular surfaces (such 
as branes) would not allow variations of the position of the surface [31, 32]. If this folklore were true then it would be 
impossible to find an action principle in a doubly covariant form. In the codimension-1 case, however, this difficulty 
was ameliorated and actually resolved by introduction of Lagrange multipliers imposing continuity of induced metric 
[16]. It should be possible and worthwhile to extend the action principle to the codimension-2 case. 
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It is useful to decompose the perturbed quantities by the scalar, vector, and tensor harmonics in Minkowski space. 
For the time being, the number of dimensions is taken to be n. rj^dx >1 dx 1 ' denotes n-dimensional Minkowski metric. 
Indices are raised or lowered by rj. 

The scalar harmonics are defined as Y — cxp(ifc M x A1 ), by which any scalar function of the position in the Minkowski 
space, f(x), can be expanded as 
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(Al) 
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where dfc = Il^=o dk^. 

The vector harmonics arc defined as 



V(tv = V> V (L) ^ = d,Y = ik,Y, (A2) 

where an arbitrary constant vector satisfies k^u^ = for non-null k^, i.e., k^k^ ^ 0. For null fc M , it satisfies 
r* 1 !/^ = for an arbitrary timelike vector t^ 1 in addition to the above condition. By the harmonics, any vector field 
v fl (x) can be expanded as 

«m = J dk ( C (T) {k)V (T)li + c (L) (fc) V (L)/1 ) , (A3) 

where the transverse part satisfies J dfcc( T )V( T ) M = 0. 
The tensor harmonics are defined as 



T(tl)hv = d,j,V(T) v + d„V(T)p = i {u^K + u u k^) Y, 

2 / 2 \ 

T(ll)vv = 9 M V( L) , y + dvV^n - -d a V( L ) a rin V = I -2fc M fc„ + -k^k^rj^ j Y, 



(A4) 



s Miy is an arbitrary transverse-traceless constant symmetric tensor that satisfies k^s^ v = s% = for non-null k^, i.e., 
fc^fc^i 7^ 0. For null k^, it satisfies t^s^ v = for a timelike vector in addition to the above condition. By the 
harmonics, any tensor field t^ v {x) can be expanded as 

V = J dk ( C (T) (k)T {T) ^ + c (TL) (fc)T (TL)My + c (LL) (fc)T (LL)AJ „ + c ( Y) (k)T {Y) ^) . (A5) 



APPENDIX B: BULK GAUGE TRANSFORMATIONS 



For our analysis it is useful to decompose the perturbations according to the transformation property from the 
viewpoint of four-dimensional space. Now all the geometrical quantities are decomposed by the harmonics on the 
Minkowski space summarized in Appendix A. The perturbations of the metric and gauge potential are, respectively, 
decomposed as 

Sg M Ndx M dx N = J d/c((/i (T) T (T>i , + h {Th) T {TL)fll , + h {LL) T {hL) ^ + h {Y) T {Y) ^) dx^dx" 

+ ( V)a V)M + fyL)aV(L) J dx a dx^ + h ab Ydx a dx b ) (Bl) 

SA M dx M = J dk ((a (T) V (T)(U + a ( L)V( L)([1 ) dx^ + a a Ydx a ) . 

Indication of integration is often omitted in the main text. 

By the double covariance shown in the main text, we can arbitrarily transform the bulk gauge. Infinitesimal 
coordinate transformation and U(l) gauge transformation for the perturbations are, respectively, denoted as 

SgMN — » 5g M N - VmCw ~ ^n£m, SA m — » SA M + d M X- (B2) 
Gauge parameters are decomposed with the harmonics as follows 

£ M dx M = [dk + l(L)V( L ) M ) dx» + t a Ydx a ) , 

J r (B3) 

X = / dfcxF. 

We assume that the bulk perturbations are axisymmetric. In such a case f's can only depend on r. Note, however, 
that the U(l) gauge parameter x still can vary with <j>. By these gauge parameters, the metric components are 
transformed in accordance with their types. 
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The tensor part is gauge-invariant 

h(T) -> h (T) = h (T) . (B4) 
As for the vector part, components are transformed as 

-(TL) - £(T)j 

(B5 ) 





■* fr(TL) 


= ft(TL) 


-|(T), 


h(T)r ~ 


-> ft,( T ) r 


= tyT)r ■ 


- £(T) + 


h( T )4, ~ 


"» V)0 


= ^(T)0- 




O(T) - 


■* 0( T ) = 


= a (T ), 





where the prime denotes derivative with respect to r. Setting £(t) — ^(tl) gives /i(tl) = 0, which is the gauge 
condition utilized in section IIIB2. 

As for the scalar part, components are transformed as 

ft(LL) — » ft(LL) = ft(LL) - £(L)j 

2 - 

^(L)r - * fyL)r = ^(L)r ~ 6(L) + ~£(L) - £r, 
^(L)0 — » ^(L)0 = ^(L)0 - 

/' " 

>■ 2^^, ~^~^,T-) 

h h -h -i< ( B6 ) 

— * h<t,4, — — // £ r , 
a (L) -» a(L) = a (L) + X - -4y , 

a r — » a r = a r + <9 r x - A > 

£(L) = ^(LL), 

2 

£r = /l(L)r + _/l (LL) - ^(LL) > 

^ f Jc f(r') ' 



Setting 



X = -0(L)+ A 



with C" an arbitrary constant gives /i(LL) = ^(L)r = = a (L) = 0. This is the gauge condition used in section IIIB 3 
together with the following redefinitions 

h(Y) = *, h{L)4> = ft(L)^j /*rr = ($1 + $2)//, ^ = - ($1 + 3$ 2 )/, a r = «r, 5^=00- (B8) 

Note that the gauge transformation considered in this appendix is for modes with fc M ^ 0. Modes with = can 
be included in the background without loss of generality. 

APPENDIX C: THE a -> 1 LIMIT 

When the four-dimensional geometry is the Minkowski, the spacetime is given by 

Ar 2 

g M Ndx M dx N = r^dx^dx" + J^ + f{rW, A M dx M = A(r)d0, (CI) 
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where 

/ = - 7! - T^6' A = ^- (C2) 

Parameters (/ib, b) are replaced by the roots (r+, r_) satisfying /(r±) = 0. In the r_/r + — > 1 limit, the metric function 
becomes 

/ = 2A 6 6 2 ^1 - (^-^j ^ + 0(r r ) 3 , (C3) 

where r = (r + + r_)/2 and e = (r + — r_)/2. Then a coordinate transformation 

w= ;r__r^ ^ = 2A 6 e0 (C4) 
gives the following representation to the metric: 

SMivda^d^ ~ r^d^dx" + -L f/»d«; 2 + , (C5) 

2A 6 V f{ w )J 

where / (to) = 1 — io 2 . The geometry of the extra space is locally 2-sphere with the radius R = 1/ V2A 6 . The gauge 
potential is transformed as 

The constant can be always removed by some U(l) gauge transformation. Then, in the e — > limit, we obtain 

A M dx M ~ RA(w)d<p, (C7) 

where = — iu. 

APPENDIX D: REGULAR QUANTITIES 

In this Appendix, we show independent regular quantities of each type of perturbations on the a = 1 background. 
See Appendix B for our gauge choices. These quantities are evaluated in the inertial frame on the brane and the 
values are required to be regular. In our position, regularity has been considered to work complementarily as boundary 
conditions. 

As for tensor-type perturbation, independent quantities in the inertial frame on the boundary are h(r) itself and, 
for example, 

(e^ M (e w ) p (e u ) N (e w ) Q SCMPNQ = -\ (wh{ T) + - m 2 ) h {T ^j T (T) ^, (Dl) 

where {e a } denotes (unperturbed) vielbein. From this combination, we deduce that 

hp), h[ T) (D2) 

must be regular at the boundaries. 

As for vector-type perturbations, independent quantities in the inertial frame on the boundaries are, after using 
equations of motion, 

(e AI ) M (e w ) p (e v ) N (e v ) Q 8C M pnq = ^_ 2 ft(T) y ^(TL)^ + ^(T)^ ^CiVi 

{e„) M {e v ) N 5F MN = a (T) {d„V (T)v - d v V (T)ll ), (°3) 
(e w ) M (e^) N SF M N = y/1 ~ w 2 a{ T) V {T)ll . 
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From this combination, we deduce that 

h 

1 — w 



h[ T)ip , o ( T), y/l -w 2 a[ T) (D4) 



must be regular at the boundaries. 

As for scalar-type perturbations, independent quantities in the inertial frame on the boundaries are, after using 
equations of motion, 

5R = -^Cu 2 $i + 4w$' 2 )Y, 

5 (R MN m>n>R M ' N ' mn) = -2[/i 2 (5$i + 6$ 2 ) + 20 W $ 2 ]r, 

1 (D5) 
(e w ) M (e^ p (e w ) N (e v f6C MPNQ = --^ 1 + ^ 2 )Y^ (p^u), 

(e,) M (e p ) p (e v ) N (e w f5C MPNQ = ~ (vT^^ + 4$^) - — ^=(<t>i + 2$ 2 )) v^p (/i, 1/ ^ p). 



From this combination, we deduce that 



2w 



$ 2 , $' 2 , $!, v 7 !-^ 2 ^- ; - ($i+2$ 2 ) (D6) 



V 7 ! 
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must be regular at the boundaries. 



APPENDIX E: ACTION FOR THE BACKGROUND BOUNDARY CONDITION 

In this Appendix we derive the action for a codimension-2 brane. We start with a general axisymmetric metric 

g MN dx M dx N = g al3 {x)dx a dx p + L(x) 2 d<p 2 , (El) 

where, throughout this Appendix, a and (3 take a, f3 = 0, . . . , D and x a does not include 4>. The axis of rotation is 
defined by L{x) = and the period of the angular coordinate cf) is A<fi. We define hm as the out-directed unit normal 
to each constant L surface 

n M = . = : = -rr-r, E2 

y/V M LV M L d±L 

where d±L = n M 8mL. Introducing radial coordinate by dp = nudx M , the (D + l)-dimensional metric g further 
admits a decomposition 

g a0 dx a dx' 3 = q pu (y, p)dy»dy" + dp 2 . (E3) 

It is clear that the above defined vector field nu is identical with that given in the main part as orthogonality to other 
vectors is reproduced. Further, since ip M riM oc dL/dcf> = and e^ M riM dL/dy^ = 0, we see that L only depends 
on p. This means that each hypersurface of constant L corresponds to that of constant p. We shift the origin of p 
without loss of generality so that p — > 0+ corresponds to L(p) — > 0+. 

Now a neighborhood of the axis of the symmetry L = is foliated by an one parameter family of (D + l)-dimensional 
hypersurfaces E^ each placed on constant p — £ > 0. We also define D-dimensional submanifolds Be contained in E^ 
by restricting 4> = const. The metric of Be is <7 M „ ($/,■£). 

We would like to place a codimension-2 brane at p — but, in order to regularize the system, first consider the 
hypersurface and define the physical brane as the I — > limit of E^. In this limit the induced metric on the brane 
is given by = q^ v {y, 0). Thus, we assume that q^ v and its curvature tensors remain regular at p = 0. 

Let us consider one such hypersurface E^. The whole manifold M is seen as union of two parts M. = M. + U M.~ 
such that M + and M.~ corresponding to p > £ and p < £, respectively. The Einstein-Hilbcrt action is decomposed 
to have two Gibbons-Hawking terms associated with two sides of T,e as 

S EH = / d D + 2 x^-gR 



M 



ME+2 



f d D+2 x^R-2 f d D+l x^qLK+ f d D+2 Xsf=g~R + 2 / d D+1 xV=qLK) . (E4) 

J M+ JdM+ J M~ JdM- ) 
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Note that in this case the direction of the normal is from J\4~ to A4 + . Trace of the extrinsic curvature of a constant 
p surface is given by 

k = -La M (V~gn M ) = %^ + = w K + k. (es) 



In the t — > 0+ limit, the Gibbons-Hawking term of (9.M + has the following limit 

/ d D+1 x^LK= [ d D y^ ( * d<f>L (^K + k) -» / d D y,/=q~A<f>d ± L, 

JdM+ JBt JO ^ ' JB 



(E6) 



where we regarded if as a finite quantity in the limit. On the other hand, the topological terms associated with M 
and dM." give 

/ d D+2 x^R + 2 / d D+1 x^LK 

JM- JdM- 

= J d D y^=q~^ ^ d$ JdpL^R + R} + 2^ d0L(«X + ^ , (E7) 

where we have decomposed the (D + 2)-dimensional Ricci scalar R into two-dimensional (p, 4>) part ( 2 >R and other 
parts R. In the £ — > 0+ limit ( 2 ^i? ~ 0(1/L), whereas -R ~ 0(1) because of the regularity of the induced metric 
Qn^(y, P = 0) on the brane. 

Therefore, by the Gauss-Bonnet theorem, terms in the square brackets give 



J *d<t>JdpL (Wr + R} +2^ *<tyL(( 1 >Jr + #) -» jj d(f>dpL^R + 2 J d(/)L^K = 4tt. 



(E8) 



Then the Einstein-Hilbert action becomes 



2 

The matter action is of the tension of the brane 



( I d D+2 x^R+ f d D y^(-2Acj ) d ± L + An)) . (E9) 
Jm Jb J 



^matter 



/ d D y^—qa. (E10) 
Jb 



Variation q^ u — > g' 1 " + f)^" to the total action S = S'eh + ^matter gives the boundary condition as an equation of 
motion 

- 2A<f>d ± L + 4tt - = 0, (Ell) 

which is actually the same as (2.11). 
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